Analytical Representation of a Black Hole Puncture Solution 
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The "moving puncture" technique has led to dramatic advancements in the numerical simulations 
of binary black holes. Hannam et.al. have recently demonstrated that, for suitable gauge conditions 
commonly employed in moving puncture simulations, the evolution of a single black hole leads to a 
well-known time-independent, maximal slicing of Schwarzschild. They construct the corresponding 
solution in isotropic coordinates numerically and demonstrate its usefulness, for example for testing 
and calibrating numerical codes that employ moving puncture techniques. In this Brief Report we 
point out that this solution can also be constructed analytically, making it even more useful as a 
test case for numerical codes. 
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Numerical relativity simulations of binary black holes 
have recently achieved a remarkable break-through. 
Since Pretorius's initial announcement of a successful 
simulation of binary black hole coalescence and merger 



ll^severalgroups have reported similar success (e.g 

1 IE IE U)- Of these, 0, @ adopt a "generalized har- 
monic" formulation of general relativity [3, [lfj and elim- 
inate the black hole singularity from the numerical grid 
with the help of "black hole excision" . All the other 
groups adopt the BSSN formulation of the ADM equa- 
tions [ll|, OjJ together with a "moving puncture" method 
to handle singularities. 

The idea of the original puncture method was to fac- 
tor out from the spatial metric (or more specifically from 
the conformal factor) an analytic term that represents 
the singular terms at a black hole singularity, and treat 
only the remaining regular terms numerically. This ap- 
proach is v ery successful for the construction of initial 
data (e.g. [l3|]), but did not achieve long-term stable 
evolutions in dynamical simulations (e.g. [3, [IB]). The 
problem may be associated with the need for a coordi- 
nate system that leaves the puncture - and hence the 
black hole singularity - at a pre-described location in the 
numerical grid, given by the singularity in the analyti- 
cal function. The break-through in the recent dynami- 
cal puncture simulations is based on the idea of using a 
"moving" puncture in which no singular term is factored 
out. Care is taken that the singularity never hits a grid- 
point in the numerical grid, but otherwise the puncture is 
allowed to move around freely. With a set of suitable co- 
ordinate conditions, found empirically, this prescription 
leads to remarkably stable evolutions. Clearly, this raises 
the question how it can be that the presence of singular- 
ities does not spoil the numerical calculation. This issue 
has been clarified recently by Hannam et. al. [1(1 [13 • 

For a single black hole, a moving puncture simulation 
starts out with a slice of constant Schwarzschild time ex- 
pressed in isotropic coordinates. These coordinates do 



not penetrate the black hole interior, and instead cover 
two copies of the black hole exterior, corresponding to 
two sheets of asymptotically flat "universes", connected 
by an Einstein-Rosen bridge at the black hole horizon. 
The singularity at isotropic radius r = 0, where the con- 
formal factor i/j = 1 + M/(2r) diverges with 1/r, cor- 
responds to the asymptotically flat end of the "other" 
universe, and is therefore a coordinate singularity only. 

Typically, moving puncture simulations use some vari- 
ant of the "1+log" slicing condition [l8| 



(dt - Pd^a 



-2aK 



(1) 



as well as a "f -freezing" condition [IS]. In [16[ the 
authors show that with these coordinate conditions, 
the evolution quickly settles down into a new time- 
independent solution that is different from the initial 
data. In this new time-independent solution the confor- 
mal factor features a 1/ yfr singularity at r = instead of 
a 1/r singularity, meaning that the slice has disconnected 
from the other asymptotically flat end, and instead ter- 
minates on a surface of finite Schwarzschild coordinate 
radius. For the slicing condition ([T]), [l6| found this lim- 
iting surface to be at a Schwarzschild (or areal) radius 
of R = 1.3M. The numerical grid therefore does not in- 
clude the spacetime singularity at R = 0. Instead, the 
singularity at r = is again only a coordinate singular- 
ity, which helps to explain the success of this numerical 
method. 

In [l7| (hereafter HHBGSO) the authors consider the 
alternative slicing condition 
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Dynamical evolutions again settle down quickly into a 
new time-independent solution, which in this case has 
to be a maximal slice (with K = 0) of Schwarzschild. 
A family of time-independent, maximal slicings of the 
Schwarzschild spacetime is given by the following expres- 
sions [20| for the spatial metric 
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the lapse 



j ij dx i dx j = f- 2 dR 2 + R 2 d 2 VL, 



a = f, 



(3a) 



(3b) 
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and the shift vector 



R 2 ' 



(3c) 



Here we assume spherical polar coordinates, and the 
function / 



/= 1 
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depends on both the areal radius R as well as an arbitrary 
parameter C. For C = we recover Schwarzschild coor- 
dinates. HHBGSO provide evidence that the dynamical 
evolution, starting from a slice of constant Schwarzschild 
time and using the slicing condition ([2]) , settles down to a 
member of the family ([3]) with C = 3\/3M 2 /4, which has 
a limiting surface at R = 3M/2. For numerical purposes 
it is desirable to obtain this solution in isotropic coor- 
dinates. HHBGSO construct this solution numerically, 
and then proceed to show that the solution is indeed 
time-independent. This solution provides a very power- 
ful test for numerical codes that use the moving puncture 
method. Using this solution as initial data, the time evo- 
lution should lead to a time- independent solution, so that 
any deviation from the initial data is a measure of the 
numerical error. 

In this Brief Report we point out that this solution 
in isotropic coordinates can be constructed analytically, 
making it an even more useful tool for code testing. 

To transform the solution ([3]) into isotropic coordi- 
nates, we identify the spatial metric (f3"a|) with its coun- 
terpart in isotropic coordinates, 



f- 2 dR 2 + R 2 d 2 n = ^{dr 2 + r 2 d 2 fl), 



(4) 



where ip is the conformal factor and r the new isotropic 
radius. From this identification we have 



R 2 = ^r 2 



and 



2 dR 2 



ip 4 dr 2 , 



which we can combine to find 
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Vi? 4 - 2MR 3 + C 2 
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In the general case, the right-hand side may be expressed 
in terms of elliptic integrals of the first and third kinds, 
which of course cannot be written in terms of elementary 
functions. The integral simplifies, however, when two or 
more roots of the quartic equation i? 4 — 2MR 3 + C 2 = 
coincide, which occurs precisely when the discriminant 



16 C* 4 (16 C* 2 - 27 M 4 



(8) 



vanishes (see, e.g., [21]). This happens for two different 
values of C 2 . In the case C = 0, integration of Q yields 
the familiar isotropic form of the Schwarzschild metric. 




FIG. 1: The conformal factor tp as a function of isotropic ra- 
dius r (solid line). The short-dashed lines shows the asymp- 
totic value tp = a/3M/ (2r) in the limit r — > 0; the long-dashed 
line shows the asymptotic value ip = 1 + M/(2r) in the limit 
r — > <x. The insert shows the same functions on a log-log 
scale. 



Remarkably, the discriminant also vanishes for C = 
3\/3Af 2 /4, the case of particular interest here. In this 
case, the quartic polynomial has a double root at R = 
3M/2, and the integral can be written 
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(R - 3M/2) ,/R 2 + RM + 3M 2 /4 ' 
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which yields 
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2R - 3M 



8R + 6M + 3(8i? 2 + 8MR + 6M 2 ) 1 / 2 
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2R + M + (4JT + AMR + 3M 2 ) 



For the positive sign, we then have 

'2R+M+ (4i? 2 + AMR + 3M 2 ) 1 / 2 
_ 1 

(4 + 3y / 2)(2 J R- 3M) 

8R + 6M + 3(8i? 2 + 8MR + 6M 2 ) 1 / 2 

M 2 



const. 
(10) 
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R 2R 2 
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where we fixed the constant of integration by requiring 
that r — > R as R — ► oo. Eq. (fTTj) gives the value of the 
isotropic radius r as a function of the areal radius R. As 
expected, we find r — > in the limiting case R — > 3M/2. 
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FIG. 2: The lapse a (top panel) and the shift /3 r (bottom 
panel) as a function of isotropic radius r (Compare Fig. 5 in 
HHBGSO). 



From eqs. ([5]) and (fTTj) , we obtain the conformal factor 



AR 



1/2 



2R + M + (AR 2 + AMR + 3M 2 ) 1 / 2 

8R + 6M + 3(8R 2 + 8MR + 6M 2 ) 1/2 
(4 + 3\/2)(2-R-3M) 



(12) 

1/2V2 



Although eq. (|11|) cannot be inverted to yield R as a 
function of r, eqs. (fTTj) and (112|l may be used parametri- 
cally to make a plot of ^ vs. r (see Fig. [1]). Figure Q] also 
includes the asymptotic limits 
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The lapse 
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and the isotropic shift 



dr 
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Rf R 2 
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may also be plotted vs. the isotropic coordinate r (see 
Fig. [2]) by expressing all quantities parametrically in 
terms of R (compare Fig. 5 in HHBGSO). 

As demonstrated in HHBGSO, this solution for the 
conformal factor rp, the lapse a and the shift (3 r is time- 
independent when evolved with the slicing condition (jlj 
for the lapse a and a T-freezing condition for the shift 
(3 l , which, in this case, leaves the shift constant (see 
HHBGSO for details). Under these conditions this so- 
lution represents the asymptotic state of a moving punc- 
ture evolution. This solution therefore provides a very 
useful test for moving-puncture codes: when adopted as 
initial data, any deviation from these initial data during 
the course of the evolution represents a numerical artifact 
and hence a measure of the numerical error. We believe 
that the analytic form of this solution presented in this 
Brief Report will significantly simplify the implementa- 
tion and evaluation of such code tests. 
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